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We calculate the leading Weyl term of the counting function for a mono-crystalline quartz sphere. 
In contrast to other studies of counting functions, the anisotropy of quartz is a crucial element in 
our investigation. Hence, we do not obtain a simple analytical form, but we carry out a numerical 
evaluation. To this end we employ the Radon transform representation of the Green's function. 
We compare our result to a previously measured unique data set of several tens of thousands of 
resonances. 
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I. INTRODUCTION 

In solid state physics quantum chaos [3- & 0- and other applications a smooth approximation to the 

level density or, equivalcntly, to its integral referred as the counting function (or staircase function) is often needed. 
Typically such an approximation is obtained by semiclassical methods |(| . Alternatively, related techniques stemming 
from geometrical quantization and index theorems Q are occasionally employed. In the context of quantum chaos, 
this problem of counting the number of states has been studied thoroughly for finite systems, in particular, billiards. 
The latter are realized in microwave resonators. Quantum chaos methods were successfully applied to elastodynamical 
systems 0, H, 0, H^, El 113 ■ The semiclassical limit corresponds to a ray limit in which the wave lengths are 

much smaller than characteristic length scales of the resonating device. However, anisotropy enters as an important 
feature of these systems which is not present in scalar Schrodingcr quantum mechanics. The problem of anisotropy 
also shows up in the calculation of the heat capacity in solid state physics which directly relates to the counting 
function. 

In this contribution, we study a mono-crystalline quartz sphere. We calculate the first approximation to the 
counting function which is referred to as the Weyl term in semiclassics. We combine analytical with numerical com- 
putations. Furthermore, we have a unique data set at our disposal which was measured previously with extraordinary 
resolution |lrll |. It comprises tens of thousands of resonances, enabling a comparison with our theoretical findings. 

Even though the particular shape of a sphere was chosen the system will not be rotationally invariant due to the 
underlying medium, a— quartz, which is anisotropic. This anisotropy leads to non-separabillity of the wave equation. 
This makes the problem so complex that there is no way of calculating spectra than by brute force numerical methods. 

A major motivation for the present investigation is the renewed interest in elastomechanical problems which are 
increasingly important for micro-electro-mechanical systems (MEMS). In particular, a better understanding of spec- 
tral properties is expected to improve the control of the resonance and scattering pole structure of MEMS. This in 
turn would allow designing better filters and actuators [l^j. Thus deeper insight into the role of anisotropy is called 
for. 

The article is organized as follows. In Sect. [HJ we compile general features of elastodynamics as needed in the 
present context. We work out the Green function in Sect. IIIll The data accumulation is briefly reviewed in Sect. II Vl 
we compare to the experimental findings in Sect.0 We summarize and conclude in Sect. IVT1 
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II. ELASTODYNAMICS 

In Sect. Ill A1 we sketch some properties of elastic waves in anisotropic materials. We discuss the dispersion relation 
and the slowness surfaces in Sect. Ill Bl The Green's function is worked out in Sect. Ill CI The Einstein summation 
convention is assumed throughout. 

A. The elastic wave equation 

In linear elastodynamics in the time domain, the wave equation reads [Tsl IT^| 

d 2 ui _ d 2 Uj 

ClJkl d Xj dx k ~ 9 dp ' [1) 

where the displacement field itj describes the increment of a fictitious point particle at its position (x\, x%, X3). The 
quantities c.^i form the elasticity tensor and p is the mass density. For our system, we may assume homogenity of 
the material, that is, the mass density and the elasticity tensor are constants. In an isotropic material, the elasticity 
tensor has only two independent entries which yield, together with the mass density, the longitudinal velocity for 
pressure waves and the transverse velocity for shear waves. In our case of anisotropic quartz, the velocity of an elastic 
wave also depends on the direction in the crystal. This means that the elasticity tensor has six independent entries. 
The elasticity tensor always fulfills the following symmetry conditions 

Cijkl = Cijlk = Cfdij ■ (2) 

The wave equation Q is valid in the regime where the linear relation 

a ij = c tjklUkl (3) 

holds between the stress tensor cr^ and the strain tensor 

1 / duk dui \ 



Equation is the most general form of Hookc's law. The wave equation in the frequency domain which corresponds 
to Eq. reads 

where the angular velocity is denoted u). 

We are interested in the elastic vibrations of a finite objects, the elastic resonator, which is confined by a surface. 
Matching the situation in the experiment we assume free boundary conditions, that is, there is no normal stress on 
the surface. This is equivalent to the condition 

= c i: jki u k i rij (6) 

with the strain and normal vector nj. The boundary condition © is of Neumann type and makes it possible 
to mode convert at the boundary. In an isotropic material, a longitudinal (transverse) wave hitting the boundary is 
reflected and generates, apart from some special situations, a second wave which is transverse (longitudinal). These 
waves have different velocities and leave the point of incident under different angles. This formally corresponds to 
quantum mechanics for a free spin-one particle whose effective mass depends on the spin degrees of freedom. The 
particle is confined in an enclosure, corresponding to the resonator. Not surprisingly, mode conversion shows a much 
richer phenomenology in an anisotropic material. 

B. Dispersion relation and slowness 

We consider a plane wave with wave vector kj , angular velocity lj and polarization vector pi , 

u[ pw \xi,x 2 ,x 3 ) = pi exp (ikjXj - iwt) . (7) 



FIG. 1: Quasi- Longitudinal slowness. 



Inserting this into the wave equation yields the dispersion relations 

Cijkikjkk pi = paj 2 Su pi . (8) 

(Here kj and arc the j th and th fc-th component, respectively, of the wave vector.) The wave numbers and the 
angular velocity occur in equal powers. This is exploited by considering their ratio, the slowness |20j| 

Si =ki/u . (9) 

The dispersion relations J5J) then motivate the introduction of the matrix 

Su = --Cijki SjS k + Su . (10) 

Thus, Eq. © can be cast into the form Supi = 0. For solutions to exist, one must have 

detS = 0. (11) 

This defines the slowness surface. As each entry in Eq. (|ll)|) is quadratic, the full determinantal condition becomes 
a polynomial of order six in the variables Sj. This condition gives rise to a surface depicted in Figs. ^ 121 an d ED 
The elastic constants c^ki for quartz have been used, see Sec. The slowness coordinates are, except for a constant 
factor, coordinates in the space dual to the configuration space of the resonator. The units on the figures are in inverse 
velocities corresponding to velocities of the order 5000m/scc. The surface consists of three sheets, S = U„ =1 E Q , 
where only the inner is convex. One clearly sees the i^-symmctry of the crystal, that is, after a rotation of 27r/6 around 
the symmetry axis a reflection in the plane orthogonal to the axis transforms the surface into itself. This surface is 
the remains of the dispersion surface after the frequency has been scaled out. Thus a fixed direction corresponds in 
general to three points on the slowness surface. The inner sheet is associated to an almost longitudinal polarization 
whereas the outer sheets are associated with the transverse polarization. Hence the terms quasi-longitudinal and 
quasi-transverse are used. 

As already pointed out, our elastodynamical problem has formally much in common with quantum mechanics for 
a particle with spin confined to an enclosure, which could be viewed as a billiard. Inside the billiard, the particle 
moves on straight lines and it mode converts upon reflections at the boundaries. This is the situation in configuration 
space. Employing slowness, we have mapped the problem of one spherical billiard for a particle with spin onto the 
problem of a particle without spin (in a dual space) simultaneously moving in three non-spherical billiards, that is, 
in the spaces confined by the three slowness surfaces. A trajectory in configuration space consists of a sequence of 
straight pieces with a certain modal character. In slowness space, the particle jumpes between the billiards which are 
specific for a given modal character. 



C. Group velocity and polarization 

The Green's function to be introduced later on is a matrix due to the vector character of the elastic field. Thus, 
it is useful to work with a proper basis in matrix space. Such a basis is formed by the projectors e a for a given 
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FIG. 2: Quasi- Transverse 1 slowness. 




FIG. 3: Quasi- Transverse 2 slowness. 



polarization a. We shall consider the case where wave vectors belong to a fixed sheet of polarization a and discuss 
the associated group velocity. Below we denote polarization indices with Greek letters whereas geometrical indices 
are in Roman letters. 

A projector e a is a 3 x 3 matrix which we introduce, corresponding to Eqs. jEJ and 1)10(1 . by writing 

k Q • c/p ■ k Q = uj 2 e a , (12) 

that is Cp qrs / pk a ^ p k ayS = u> 2 e a . qr . Each projector corresponds to one of the three possible polarizations a. The 
projectors e a satisfy the completeness relation 

^e Q = l 3x3 (13) 

a 

and the relation 

e Q • ep = 8 a p e a , (14) 
where the dot indicates matrix multiplication. The dispersion relation for a fixed polarization becomes 

(k Q • c/p ■ k Q ) : e Q - uj 2 = 0, (15) 
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FIG. 4: Twisting of polarization: Shear 1 



where the contraction of all indices is denoted by the :-sign. 

A natural quantity that now enters the discussion is the group velocity for a given polarization a 



= dk~ a ' (16) 



describing the change of angular velocity with the wave vector for that polarization. Varying Eq. i|15|) with respect 
to the wave vectors and the angular velocity 

(dk Q • c/p ■ k Q ) : e Q — du u = (17) 

and finally scaling out the angular velocity gives 

du . . . . . 

v g , tt = ^r— = (c/p ■ s a ) : e a . (18) 

Above, also the projectors could have a variation which we have not taken into account. However, since 

(k Q ■ c/p ■ k a ) : de a = u 2 e a : de a (19) 

and 



2 e Q : de Q = d(tr (e Q ■ e Q )) = d(tr (e Q )) = d(l) = , (20) 

this variation does not contribute. 

From Eqs. © and (|15fl we see that s Q and v gjQ , are mutually polar reciprocal to each other, that is, 

1 = (sq • c/p ■ s a ) : e a = s a ■ (c/p ■ s a ) : e a = s a ■ v g , Q . (21) 

This relation can be used to construct the group velocity and holds in general for systems which are homogeneous in 
ki and lu. This condition defines the individual sheet, and in particular the group velocity is normal to this sheet. 

However, interesting complications arise as the degeneracy of this matrix is not constant. Hence there exist directions 
for which the phase velocity of the two intersecting sheets is the same. Note that this degeneracy occurs in slowness 
space. Interestingly such mode conversion points may also occur purely in x-space, e.g. in Landau-Zener problems or 
in nuclear/chemical reaction physics. In general, however, there are three distinct polarizations, and from Eq. i|lU|) 
these only depend on the direction s Q = s Q /,s Q . Furthermore for the opposite direction the same polarization is 
found. We have depicted one of the quasi-transverse polarizations on Figs. 0] and The behavior of the two 
shear components are similar with a twisting near the poles. The quasi-longitudinal, however, is much simpler with a 
polarization vector which looks like a smooth deformation of the field of normal vectors on the sphere (a hedgehog), 
see Fig. |H] 
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III. THE GREEN'S FUNCTION 



We discuss a result based on the Radon transform. Numerically it is convenient to express the Green's function as 
a certain integral over the unit sphere. For a theoretical interpretation this integral is transformed to a similar one 
over the slowness surface. Using the Radon transform the result in the frequency domain is |2lLl22|| 

G + (x,w) = G^(x,co) + G s (x,w) (22) 

with superscripts R and S indicating regular and singular part and subscript + for causal. Here we have changed 
dimensions compared to Refs. |2ll l22l | by solving for the Green's function for Eq. JSJ) where p has been divided out. 
For convenience a derivation is discussed in App.^ The regular part is an integral over plane waves: 

G« M = ^[ £^ e *.l~ldn. (23) 

Here k a ,c a and e Q refer to wave number, phase velocity respective polarization projector for polarization number a 
for the plane wave in question. The singular part 

G s (x,lu) = / r -1 (n) 5(n • x)dO (24) 

is purely real and corresponds to the static part of the Green's function. The matrix T is defined by 

Ciklj 



Fij = n k ni (25) 
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for which 

r(n)=E^°- ( 26 ) 

a 

We notice that Q25|l follows from l|15|l with k = nk and c = uj/k. 

The integrals above can be formulated in a more intrinsic form. Thus instead of integrating over the unit sphere 
S 2 the dispersion surface can be applied. Using slowness, Eq. I|23|) becomes, 

G*(x,c) = -\c/ ££^ e M--*idn. (27 ) 

Since the wave number and the angular velocity occur homogeneously in the dispersion relation the corresponding 
slowness surface is used. Applying polar reciprocity l|21(l . 1 = s • v s and denoting 9 the angle between s a and v gjC[ , 
the element of area do~ a on this surface is 

cos 9 , da a 

dQ = -^-da a = 21_ (28) 

so 

G ?( x ^) = i-E/ 2^-e*-l-*l*r a . (29) 
This can be written in a causal form as an integral over forward pointing wave vectors 

G*(x,o;) = ^wW — e*"— d<7 a . (30) 

The general result including the anti-causal Green's function reads 



G£(z,y) = ±iu,-^W da a -^e^ sr , (31) 
( 27r ) „ Jse-z± 2v s,a 

where T,^ is the positive or negative of the surface S Q , respectively. For the causal Green's function this is to be 
understood over the upper hemisphere (sn > 0). For the anti-causal Green's function, the lower hemisphere, is the 
domain of integration. Furthermore the group velocity is opposite of the direction of observation, leading to a sign 
change of the integrand. In conclusion, the discontinuity of the Green's function 

AG(i 1 y)=G + (i,y)-G.(i 1 y)=j(j ? ^r / da a e -*C*-y) (32 ) 

( 27r ) „ Jsex a 2v g,a 

is expressed as an integral over the whole slowness surface. Note here the cancellation of the (real) singular part in 
the discontinuity of the Green's function. 



IV. ACCUMULATION OF DATA 

In the experiments 0|, several shapes were studied. This article concerns the spectrum of a particular quartz 
sphere. We shall briefly sketch the experiments, further details will be given in Ref. 16]. The quartz sphere is 
supported by three transducers, see Fig. [3 The transducers excite the sphere via piezo-electric coupling. The data 
obtained in the experiment is an arbitrary constant times the complex ratio of the outgoing to the incoming electric 
signal. As the data are complex the actual phase of the response is available. The setup and the measurement 
techniques are extensions of previous work [lf| . The resolution of the setup, characterized by the Q-value, is very 
high. A typical resonance measured in experiment is shown in Fig. |SJ Due to the high Q- value ringing is observed 
because the transducer frequency changes at a constant rate. The actual position of the single resonance is to the left 
of the observed peaks as the frequency is swept upwards. Since occasionally two or perhaps even more resonances are 
close together, the observed signal appears as superpositions of such. The resonances are fitted by using a recognition 
code |23| which is based on maximum entropy principles. The code fits for each resonance the position in frequency, 
the decay width and its (complex) amplitude. It is capable of treating individual resonances with ringing and of 
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FIG. 7: Quartz sphere with supporting transducers. 

superimposing them if they are overlapping. Good performance is observed for the real and imaginary part of the 
response. 

The quartz sphere was manufactured with a high precision. The spherical shape ensures the discrete symmetry of 
the crystal, which is described by the symmetry group D3. However, D3 is broken to the symmetry group C3 due 
to the positions of the transducers in the experiment. Moreover, despite the high precision manufacturing, the shape 
cannot be perfectly spherical. Its diameter has a relative precision of 5 • 10~ 5 . At high frequencies this should lead 
to observable splittings. Finally, gravity leads to an oblate deformation of the sphere. These effects further break the 
symmetry. Hence, the vast majority, if not all, of the symmetry degeneracies are expected to be broken. 



V. WEYL CONTRIBUTION 



We now discuss the counting function. We work out the leading approximation in Sect. IV Al and evaluate it 
numerically in Sect. IV 51 In Sect. IV 51 we compare to data obtained in the measurements mentioned in the previous 
Sect.llV] 



A. Approximation to the counting function 

The smooth or Weyl part of the level density is the easiest to access. The fluctuations play an important role as 
well. In fact, the oscillating part of the level density becomes of the same order as the smooth part. Here, however, 
the focus is on the smooth part. The level density is written formally as 



w — Wi). (33) 

i 

This can be found from the discontinuity of the Green's function as: 

d(u a ) = — TrAGfa; 2 ) (34) 
2ni 

since we here use a Green's function with respect to the squared angular velocity. Now the Green's function is 
approximated with the corresponding free function. This step is discussed in |24| and below in this section and 
corresponds to that in the interior parts of the resonator the propagation is roughly free. This gives 

«^_y.„^£ J f^ (fe „_L_. (35) 

The volume of the resonator is denoted V x . The counting function 

N(u) = N(u 2 ) = / d{w' 2 ) duo' 2 (36) 
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FIG. 8: Resonance with ringing due to finite sweep velocity and high Q-value 



becomes 



3 ( 27r ) „ is£E„ v S,a 



(2tt) 3 — JseSo Vgi , 

From Eq. 1)2 1 ft and geometrical considerations of cones the volume element on a given slowness sheet is 



dV a 



1 da a 

3 V g,Q 



and therefore 



(2tt) 3 



in terms of the volumes of the slowness sheets, V a . Somewhat simpler, using frequency / 

N(f) = V x J2 V «-f 3 - 



(37) 



(38) 



(39) 



(40) 



Thus we have a re-derivation of the well-known result: namely, the number of states equals to leading order the 
available phase space volume. 

Although we deal with the anisotropic case, it is instructive to see that the formulae reduce to the result for the 
isotropic case, 



47T 

T 



(41) 



where cl and ct are the longitudinal and the transverse sound velocities, respectively. 
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TABLE I: Elastic constants for quartz at room temperature. 



units 


Cll 


C12 


C13 


C14 


C33 


C44 


lO^N/m 2 


0.868 


0.0704 


0.1191 


-0.1804 


1.0575 


0.5820 



In solid state physics the procedure is different. The resonators are unit cells in crystal lattices. This implies a 
periodicity assumption of the crystal in contrast to our case of a finite system. Lattice vibrations are connected to the 
concept of heat and the counting function for vibrations leads to the famous Debye law for the heat capacity. This is 
derived in Ref. 0. 

We will now discuss heuristically why the above derivation is nothing but the leading approximation and correponds 
to the smooth part of the counting function. Since the trace of the Green's function is taken, it is mainly the local 
propagation between nearby points which is probed. However, also the possibility of propagation via reflections from 
the boundary remains. Contributions here come from periodic orbits p but are typically of lesser magnitude than the 
direct zero length contribution. Ultimately such terms lead to oscillatory contributions to the spectral density. Thus 
expanding around periodic orbits using the method of stationary phase leads to terms of the form of amplitudes A p 
times phases 

A p exp(iwT p + ...), (42) 

where the (dimensionless) classical action k L p = uT p in terms of the length/period of the orbit controls the phase of 
the fluctuations. Averaged over a sufhcently large frequency interval the fluctuations dissappcar and what remains is 
the part from zero length orbits. The latter do not fluctuate and represents the average contribution. 



B. Numerical calculation and prediction 

At room temperature T r = 22°C, the mass density is p = 2.6485 g/cm 3 , and the diameter of the sphere measures 
87.75 mm. The elastic constants for quartz at room temperature are listed in Tab.[H They are given in Voigt notation. 
Pairs of Cartesian indices ij with i, j £ {x, y, z} are numbered as follows xx = 1, yy = 2, zz = 3, yz = 4, zx — 5 and 
xy = 6. Thus for example c xxyz is denoted C14. 

We now calculate the leading Weyl term l|40|) . The volume in slowness space is conveniently expressed as an integral 
over the unit sphere, 

V s = Y,V a = Y,\l 4^n. (43) 



n£S 2 



The three lengths of the slowness vectors are found from Eqs. (|10llll|) . In particular, the squares of the velocities are 
found from Eqs. 1(23) and 1(25^ . Therefore Eq. 1031) reduces to 

V s = \( tr(r(n)- 3 / 2 W (44) 

which is easily calculated from the eigenvalues of r(n). Finally the symmetry can be exploited: In the case of quartz 
it is enough to integrate over a fundamental region < <f> < 2tt/3 and < 9 < tt/2 and multiply by six. Using the 
elastic constants in Tab. ^ we find the volume in slowness space by numerical integration, 

V s = 1.4909 • 10~ 10 (sec/m) 3 . (45) 

This volume may be used for various shapes of resonators upon multiplication with the particular volumes in physical 
x-space. Hence, combining everything, the number of states for our sphere should grow as 

N(f)=Af + 0(f 2 ) (46) 

with the constant 

A = 5.274 • 10~ 14 sec 3 (47) 



as the result of our computation. 
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FIG. 9: Spectral staircase: experiment and fit 

The elastic constants change with temperature. In App. [BJ they are given in terms of an expansion. As the 
experiment was performed at T = 30°C, we also calculate the volumes in slowness space with these corrected constants. 
However, it turns out that the relative difference to the calculation for room temperature is of the order of 10~ 8 , 
that is, within the numerical errors. Similarly, the temperature dependence of the mass density, see App.[Bj has no 
influence cither. 



C. Comparing with experiment 

The spectrum used here was measured from 410 kHz to 960 kHz. The resonances were identified by means of the 
recognition routine [2j| mentioned in Sec. II VI Roughly 52000 resonances were found. Here the bulk term in Eq. i|46|) 
contributes with at least 43000 resonances. The staircase for the experiment is fitted with a cubic polynomial, see 
Fig. ED For the leading coefficient we find 

A fit = 5.030- 10" 14 sec 3 . (48) 

This is in good agreement with the theoretical result (|47|l . The relative error is 4.6%. 

The difference between the fitted and the experimental counting function is plotted in Fig. 1101 The sizable deviation 
between fit and experiment is striking. It resembles a sawtooth with very rapid fluctuations on top. We Fourier 
transform this oscillating structure and obtain the period spectrum, see Fig. 1111 From this we estimate the time 
corresponding to the largest peak to approximately 2 /zsec. Also other peaks at 4, 9, 13, 17, 20, . . ./xsec are observed. 
As the frequency interval sampled is approximately 500 KHz the resolution is roughly 2 /isec. It is tempting to 
interpret the largest peak as due to the shortest bouncing-ball-typc-of mode, that is, a wave packet bouncing back 
and forth along some diameter of the sphere. If we take the velocities to lie between 5 km/sec to 10 km/sec the 
shortest period is around 17.8 [isec. Therefore the peak corresponding to the sawtooth oscillations is not likely to be 
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FIG. 11: Period Spectrum 

due to such a mode and remains unexplained so far. Likewise the remaining short times less than the estimate 17.8 
/xsec are not explained. 
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VI. SUMMARY AND CONCLUSION 



We have discussed the experimental spectrum of a quartz sphere from the point of view of the spectral counting 
function. The classical phase space volume gives indeed the leading behavior of this object. On calculating this 
numerically we find it in good agreement with the experimentally fitted term. This paves the way for further studies 
of elastic spectra, such as boundary corrections to the counting function [2M |26[ and periodic orbit theory for the 
fluctuations in the level density [T3L Ibi, IT5I ] . 

Due to the boundary further resonances should be present. Less is not expected as the boundary conditions allow 
for surface waves in analogy with Neumann conditions for the Helmholtz case. It is currently an open problem to 
incorporate this surface contribution in the anisotropic case. Nevertheless such surface modes couple to the bulk and 
are not exponentially damped in this elastic case. Therefore experiments should access these states and thereby the 
surface correction to the leading volume term as well. 

The fluctuating part would involve recent semiclassical methods. Here individual closed rays, periodic orbits, control 
the oscillations in the level density. In particular this has been shown not just for quantum mechanical spectra but 
also more recent in two dimensional elastodynamic resonators. The step to three dimensions plus anisotropy has yet 
to be made. 
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APPENDIX A: RADON TRANSFORM AND THE GREEN'S FUNCTION 



Following Ref. |2l( we first solve the Radon transformed time problem. Next the solution is transformed back and 
finally Fourier transforming gives the frequency dependent Green's function. 



1. Radon transform 



The equation s = n • x describes a plane with the normal vector n in the space of the position vector x. The plane 
cuts the direction defined by n in the height s. The Radon transform /(s,n) of a function /(x) is the average 



Tlf(s, n) = f(s, n) = / f(x)S(s - n • x) dx (Al) 

over the plane for a given height s. The Radon transform of the spatial derivatives of the function /(x) fulfill 

n(Jpj(s,n)= ni (d s nf)(s,n) . (A2) 

We now apply this to the time dependent elastic problem We denote the Green's function G(x, t) and its Radon 
transform G(s,n). We find 

(r(nR 2 -Id?) G = -l S(s)S(t), (A3) 

where the matrix T is defined in Eqs. i|25[l and (|26fl . We observe that G(s,n) also depends on time t. The initial 
condition is taken as causal such that G(s, n) = for t < 0. Necessary and sufficient conditions for solving this 
problem are found by projecting onto the subspaces defined by the individual polarizations: 

(c 2 a d 2 s - d 2 t ) e a ■ G(s, n) = -e° S(s)S(t) , (A4) 

where c a is the phase velocity associated with polarization a and e a is the corresponding projector. Invoking 
d'Alcmbert's solution 

e a • G(s, n) = e Q {H{s + c a t) - H{s - c a t)) (A5) 
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with H the Hcavisidc unit step function, we find 



G(s, n) = ]T e a (H(s + c a t) - H(s - c a t)) , (A6) 
zc a 



which is the full solution. 



2. Inverse Radon transform 

The inverse Radon transform is given by 



«5(x)=-J- / dO<5"(n-x). (A8) 
8*" Jnes^ 



following from the identity 



The latter is proved using 

2tt|x| = / dCl |n-x| 

and the identities \x\" = 2(5(x), = — 47T(J(x) and A|x| = ^ . Hence, we find the Green's function in configuration 
space, 

G(x, t) = -ffi / dn V 9 t (<5(c Q t + n • x) + «5(c a t - n • x)) 

= -ITT I dn Y,^T d t S(c a t - n • x) , (A9) 

by inverting Eq. i|A6() . 

3. Frequency domain 

The frequency dependent Green's function G(x,a>) is obtained by fourier transforming (|A9() with respect to time: 
G(x,w) = 7/ d^V— / dte luJt H(t)d t S{c a t-n-x). (A10) 

The inner integral is found by partial integration: the time derivative will act either on the exponential or the step 
function leading to the regular respective singular contribution 



and 



G 5 (x, w ) = — / r 1 (n) S(n • x)dfi . 



The anti-causal Green's function is found similarly by reversing the time in say Eq. i|A6(l . In particular the same 
singular part is otained whereas the regular part is the complex conjugate of the causal regular part. 
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TABLE II: Coefficients of temperature expansion for eiastic constants of quartz. 



order 


units 


Cll 


C12 


C13 


C14 


C33 


C44 





10 11 N/m 2 


0.868 


0.0704 


0.1191 


-0.1804 


1.0575 


0.5820 


1 


xW^/K 


-48.5 


-3000 


-550 


101 


-160 


-177 


2 


xlO~ 9 /K 2 


-107 


-3050 


-1150 


-48 


-275 


-216 


3 


xlQ- 12 /K 3 


-70 


-1260 


-750 


-590 


-250 


-216 



TABLE III: Coefficients for the temperature expansion of the mass density of quartz. 



order 


units 


P 





10 3 kg/m 3 


2.6485 


1 


xW~ 6 /K 


-34.80 


2 


xW- g /K 2 


-30.04 


3 


xW- 12 /K 3 


49.08 



APPENDIX B: TEMPERATURE DEPENDENCE 



The resonator is kept at fixed temperature 

T = 30°C (Bl) 

with a precision on the order of a milli-Kclvin. This is slightly higher than the room temperature T r = 22° C, at 
which the quartz sphere was manufactured. This causes a change in the elastic constants described below in terms of 
a power series expansion around room temperature, 

Cijki{T) = Cijki,o + Cijkl,i(T - T r ) + c ijk i t2 {T - T r ) 2 + c ljkl ^{T - T r ) 3 + . . . (B2) 

We took the coefficients of this expansion from Ref. |2j], they are listed in Tab. ITU Furthermore a change in mass 
density takes place. In Tab. IIIII we list the coefficients of an expansion around room temperature, Finally, there is 
also a shape deformation, expressed in terms of the tensor of thermal expansion according to u! i = onjUj. 

The effect of temperature change is therefore mainly to change the shape and volume of the /c-part of the available 
phase space volume, namely the slowness surface. Thus the first table is most important for our purposes. We have 
calculated the new volume in slowness space at the corresponding experimental temperature to cubic order in the 
temperature change. The other, relative changes turn out to be negligible, namely of the order 10^ 8 . 
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